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Using a coarse-grained bead-spring model of flexible polymer chains, the structure of a polymeric 
nanodroplet adsorbed on a chemically decorated flat wall is investigated by means of Molecular 
Dynamics simulation. We consider sessile drops on a lyophilic (attractive for the monomers) region 
of circular shape with radius Rn while the remaining part of the substrate is lyophobic. The variation 
of the droplet shape, including its contact angle, with Rd is studied, and the density profiles across 
these droplets also are obtained. 

In addition, the interaction of droplets adsorbed on two walls forming a slit pore with two lyophilic 
circular regions just opposite of one another is investigated, paying attention to the formation of 
a liquid bridge between both walls. A central result of our study is the measurement of the force 
between the two substrate walls at varying wall separation as well as the kinetics of droplet merging. 
Our results are compared to various phenomenological theories developed for liquid droplets of 
mesoscopic rather than nanoscopic size. 



1. Introduction 

In the context of the emerging technology of microscale processing the use of surface substrates with 
micropatterns has become increasingly important for microfluidics applications. Chemically struc- 
tured surfaces that exhibit lateral patterns of varying wettability can be produced by techniques such 
as photolithography 1 ' 2 , microcontact printing 3 5 , vapor deposition through grids 6 , domain forma- 
tion in Langmuir-Blodgett monolayers 7 ' 8 , electrophoretic colloid assembly 9 , lithography with colloid 
monolayers 10 , microphase separation in diblock copolymer films 11 , etc. For patterned surfaces in 
the micrometer range, on which liquid droplets or thin liquid films are adsorbed, fascinating wetting 
morphologies have been predicted 12 19 (including "morphological wetting transitions" 12 ' 16,17 ) and 
observed 6 . Liquid bridges between droplets at walls are also of interest in the context of forces be- 
tween droplets (or bubbles, respectively) 20 ' 21 , long range forces between colloidal particles 22 ' 23 , etc. 
The theoretical treatments mostly apply phenomenological quasi-macroscopic concepts (in terms 
of interfacial tensions, contact angles, 24 ' 25 etc.); but for droplets on the micrometer scale the lack 
of knowledge on the line tension 25 36 is a serious drawback already, and the understanding of the 
density distribution of the droplet near the contact line is a difficult problem 34 ' 35 . The present drive 
towards nanoscale technology creates also more interest in droplets on the nanometer scale, and 
in fact some of the techniques mentioned above are well suited to create surface patterns on the 
nanoscale 10 ' 11 . Although sometimes macroscopic concepts do allow reasonable predictions down to 
the nanoscale 37 40 , there is no guarantee that the phenomenological theories are quantitatively accu- 
rate for bridge formation between polymer nanodroplets. We study this problem here by Molecular- 
Dynamics simulation, to complement existing knowledge on the problem by insight on a molecular 
level. 

In the next section we shall present our model, which has been used successfully in previous work 
to study static and dynamic properties of polymers in the bulk and at surfaces 37 39,41 47 . In section 
3 we shall present our simulation results on droplets on single walls, as well as on the structure of 
liquid bridges in slit pores of varying width. Section 4 is devoted to a discussion of dynamic aspects 
of bridge formation, while Section 5 summarizes some conclusions. 



2. Some comments on the model and the simulation technique 

We employ a coarse-grained off-lattice model of polymer chains, where each chain consists of N = 32 
effective monomers, which are connected by anharmonic springs. These effective bonds are thought 
to represent groups of a few successive chemical monomers along the chain, and therefore inclusion of 
torsional potentials and even bond-bending potentials is not considered. The springs are described 
by the finitely extensible nonlinear elastic (FENE) potential 

u FENE (£) = -|i? 2 Mi - ( ^~f o)2 ] , (1) 

I being the bond length, which can vary in between ^ mm and <max, and thus the equilibrium value 
£o for which Ufene{£o) = 0, and R — Imax. — io — £o ~ Anhv The s P r i n S constant K is taken as 



in previous work K/ksT = 40, and we again choose fmax = 1 as our unit of length, with R = 0.3 
(hence £ = 0.7, £ min = 0.4) 41 " 47 . 

Between the effective monomers a Morse potential acts, 

U M (r) = e M exp[-2a(r - r min )] - 2 exp[-a(r - r min )] , (2) 

r being the distance between the beads, and the parameters are chosen as r mm = 0.8, 6m = 1 and 
a = 24. Owing to the large value of a, Um(t) decays to zero very rapidly for r > r mm , and is 
completely negligible for r > 1. This choice of parameters is useful, particularly for Monte Carlo 
simulations, since it allows the use of a very efficient link-cell algorithm 41 . The Theta-temperature 
for this model is 43 feO w 0.62. We hence choose in the following a temperature fcsT = 0.49, where 
at zero pressure the system already is in the state of a dense melt, and to a very good approximation 
the gas density is zero for N = 32. Thus from a polymer droplet at a surface no chains evaporate 
into the surrounding gas. 

The adsorbing walls are treated as perfectly flat and structureless; in particular, no atomic cor- 
rugation of the walls is considered. The interactions between the effective monomers and the walls 
are represented by a Lennard- Jones potential, integrated over a (semi-infinite) substrate 39 

U w&n ( Z ) = 4^(^f - i(^il) 3 ] , (3) 

where we choose parameters cr wa y — 1, fw = 0.20 in the lyophilic part of the substrate, while in 
the lyophobic area ew = 0.05 is chosen. Typically, the lyophilic area is a circle of radius Rd, with 3 
< Rd < 15. The total lateral linear dimensions of the simulation box is chosen 64 x 64, so that (for a 
total number of N = 128 chains in the polymer droplet) it never can happen that a droplet interacts 
with its images generated by the periodic boundary condition. For a study of single droplets, the 
linear dimensions in the perpendicular direction is L = 32, and the top wall of the box is taken 
purely repulsive (omitting the attractive term from Eq.(3), and choosing also e w = 0.20). For our 
study of slit pores, we choose 5 < L < 22, and the potentials from both walls are chosen exactly of 
the same type. 

The initial preparation of droplets in equilibrium is done by Monte Carlo methods, applying the 
Velocity- Verlet algorithm and keeping the temperature constant by the Nose-Hoover thermostat 48 . 
In the case of bridging droplets, however, we resort to a Langevin thermostat which provides 
greater stability of the algorithm for small separation between the solid planes. For more details on 
our Molecular Dynamics algorithm, we refer the reader to 39 . Here we only note that due to the steep 
variation of the Morse potential a very small integration time step needs to be used, St — 0.0009 
MD time units, and typically runs over 1.1 million time steps were carried out. 

It is certainly useful to translate the time units [t.u.], used in our simulation, into seconds by 
mapping our model results to laboratory data. A typical substance used in numerous experiments 
on wetting is, for example, the PDMS (polydimethylsiloxane) melt {C2H(,OSi) n - In ellipsometric 
experiments on droplet spreading 50 one has measured a diffusion coefficient ~ 3.3 x 10~ 10 m 2 /s 

for chain lengths N = 10, and Vdiff ~ 0.3 x 10~ 10 m 2 /s for N — 20. These data can be compared 
to our measurements of precursor diffusion 39 in spreading droplets which yield ©X// ~ 0.11£o/t.u. 
in a melt of chains with N = 32. Bearing in mind that the bond length in PDMS is 51 « 1.59A, 
and assuming that, say, three chemical units form a persistent length of ~ 4A, one obtains for the 
simulation time unit lt.u. ~ 0.5ns as a rough estimate. 



3. Polymer nanodroplets on chemically structured walls: static properties 

As discussed in our previous work on droplets adsorbed on flat walls without any chemical structuring 
on the substrate, the average density distribution p(R, Z) is a function of two coordinates, the 
distance Z from the substrate (note that we orient the z-axis perpendicular to the surface and 
the origin of the coordinate system is the projection of the center of mass of the droplet on the 
substrate plane z = 0), and the distance R from the z-axis. While individual configurations of the 
droplet due to statistical fluctuations depend on the angle ip relative to the x-axis as well 37 , the 
average density distribution must have rotation symmetry around the z-axis 37-39 . When we create 
a chemical structure on the surface such that the constant ew describing the wall potential U wg j[(Z) 
in Eq. (3) has one value e w = 0.2 for < R < Rd and a smaller value e w = 0.05 outside this circle 
whereby the symmetry axis of the droplet coincides with the axis perpendicular to the midpoint of 
this lyophilic circle, of course. Therefore it is meaningful to record density profiles p(R, Z) of the 
same type to characterize the average droplet density profile as in our previous work. Fig. 1 shows 
a few representative examples. While for Rd = 15 the droplet sits fully inside the lyophilic circle 
and hence hardly differs from a droplet on an infinitely extended lyophilic surface, for Rd < 12 
the droplet always extends over the full range of the lyophilic region. One thus can see that the 
shape of the adsorbed droplet changes when the radius Rd of the lyophilic domain decreases. As 
expected, the contact angle then varies continuously with Rd, due to the interaction of the contact 
line and the boundary between the lyophilic and lyophobic regions at the substrate, and is no longer 



identical to the contact angle that applies for an infinitely extended flat lyophilic substrate (this 
contact angle is determined, for very large droplets, in terms of the polymer-wall and polymer-gas 
interfacial energies, through the Young equation 24,25 while for not so large droplets also a correction 
due to the line tension 35 ' 36 needs to be taken into account 37 ' 52 ). As discussed in our earlier work 37 , 
for nanodroplets considered here, some ambiguity in the definition of the contact-angle for these 
nanodroplets is inevitable. We follow the earlier work 37-39 , fitting a straight line to the density 
contour p(Z, R) = f in the regime 2 < Z < 4, disregarding the slight curvature of this contour 
in that region. Fig. 2 presents a plot of the resulting variation, showing that cos{9) ~ const for 
Rd > 12, while cos(8) < for Rd < 9. Comparing the density profiles shown in the insert of 
Fig. 2 with those of Fig. f one sees that indeed for Rd > 12 the droplet shape in Fig. 1 is the 
same as that for a homogeneous surface with e w = 0.2. Qualitatively, the behavior seen in Figs, f ,2 
nicely corresponds to the theoretical predictions of Lipowsky et al. 12-18 . When Rd gets very small, 
the droplet shape gradually approaches the shape that a droplet takes on a uniformly lyophobic 
substrate surfaces, as the comparison of the droplet profile for Rd = 3 in Fig. f and the droplet 
profile for e w — 0.05 in Fig. 2 shows. 

We next study the behavior of a slit pore of width L, where both walls exhibit a lyophilic domain 
of radius Rd exactly opposite to each other. When L is large enough, in equilibrium (for a fixed total 
number TV of chains) we expect that droplets containing Af/2 chains each (assuming J\f is an even 
integer) will be adsorbed exactly opposite to each other. However, when the distance L between the 
plates gets smaller, the two droplets start to interact, and a liquid bridge between both walls can 
form. This formation of liquid bridges by variation of L is illustrated in Fig. 3, choosing Rd = 14 and 
J\f/2 = 128, so that the single non-bridging droplets on the separated walls are exactly equivalent 
to the situation considered in Figs. 1, 2. One can see that for L — 22 the droplets are still separated 
and identical in shape to those seen in Fig. 1 for 12 < Rd < 15. For L < 19, however, bridge 
formation has occurred. One can clearly see the change of the bridge morphology with decreasing 
distance L between the two walls. When the separated droplets touch each other, they form an "in- 
bridge" (i.e., the curvature of the bridge surface is concave) which is typical for liquids wetting the 
substrate. This catenoid shape is displayed for the profiles when the pore width is L — 17, 15, or 13, 
respectively. On further decrease of L, the shape of the liquid bridge becomes perfectly cylindrical 
(the shapes for L = 10 and L = 9 are close to this shape), while for still smaller slit pore widths 
L (such as L = 6 and L — 5) the shape is that of an "out bridge", with a convex curvature of the 
bridge surface, i.e. a barrel-like shape. This shape is typical for lyophobic substrates, as experienced 
here by those monomers of the polymer droplet whose coordinates R(X, Y) exceed the radius Rd 
of the hydrophilic domain. Again the behavior in Fig. 3 is very nicely consistent with the behavior 
as predicted by the theory for quasi-macroscopic droplets 12-19 , and hence we again find that these 
phenomenological concepts due to Lipowsky et al. 12-18 work qualitatively down to the nanoscopic 
scale. 

In view of the fact that there is some ambiguity in the precise numerical estimation of the contact 
angle of nanodroplets, as mentioned above, it is of significant interest to estimate additional proper- 
ties quantitatively, which could be compared to analytical theories on this problem. Such properties 
are the base radius R\g^ of the droplet (Fig.4) or liquid bridge (Fig. 5), the height of the droplet 
H (Fig.4), and the midheight radius Ry & ^(Z = H/2) of the droplet (Fig. 4) or of a liquid bridge 
Ry & ^(Z = L/2), (Fig. 5), respectively. Figs. 4,5 demonstrate that these quantities can be measured 
with relatively small statistical errors in our simulations. From such data it also is evident that a 
pronounced change in the behavior occurs at Rd ~ 8.5 (Fig. 4). For a macroscopic sessile droplet 
(satisfying the Young equation with the contact angle 9 and having a sphere-cap shape) we simply 
would have the relations in terms of the sphere radius r 

r sin(0), H = r(l - cos(6»)) , 

and hence 

= (1 - cos{6))/ sin (0) » 6>/2 , 9 -f . 

Indeed this relation is roughly fulfilled for our model. 

It also is of interest to calculate the pressure in the liquid bridge, using the virial formula 46 ' 53 . 
Fig.( 6) shows that the pressure is positive for the smallest distances between the plates (L < 6) 
only, while for larger distances the pressure is negative. This observation already indicates that 
this situation is unfavorable, if we would not enforce the distance L between the plates as given 
parameter, such distances would not occur if the walls could freely move against each other. This 
fact is very clearly borne out when we compute the normal force between the walls (Fig. 7). Of 
course, for large distances L for which two separate droplets occur on each wall that do not yet 
touch the force is zero, while for bridging droplets an attractive force arises, and only for very 
small L, where the bridge is squeezed into the lyophobic part, the force becomes repulsive (Fig. 7). 
The force goes through zero at L « 6. Since it has been shown 19 that for the underlying catenoid 
geometry a closed analytical expression for the surface area dependence on the wall separation L does 
not exist, wc tentatively use a relation proposed by Swain and Lipowsky 14 for the bridge between a 
single pair of opposing lyophilic stripes. The corresponding contact angle 9* can be calculated then 



from 9* = -tan' 1 (2R D / L) w 100° and agrees with our observations. As expected, this contact 
angle exceeds n/2. The strongest attractive force actually does occur near L w 10, where the angle 
is 7r/2. For L > 12 an almost linear decrease of the force sets in, before it discontinuously jumps to 
zero for L = 22. 

It is interesting to note that very long range attractive forces have been experimentally detected 
in AFM measurements of forces between a polystyrene sphere and liquid interfaces 22,23 , showing a 
linear variation with distance on the scale of 30 nm. Unlike such experiments, in our simulation an 
electrostatic origin of such long range forces is excluded by construction of our model. The force 
seen in in Fig. 7 is entirely due to the interplay of the various interfacial interactions that control 
the shape of the liquid bridge. One should also point out that both the course of the bridging 
force against inter-plate distance L as well as the particular bridge configurations corresponding to 
various parts of the force distance curve closely match some recent results 19 on classical structureless 
droplets obtained by two different surface minimization techniques. 



4. Dynamic aspects of bridge formation 

For L — 21.5 one can observe initial states where the two substrates still carry separate droplets, 
but the wings of their density distributions already overlap, and this interaction between the two 
droplets starts a merging process which leads to the formation of a liquid bridge. Figs.( 8) and (9) 
analyze the dynamical aspects of this merging process of the two droplets in more detail. One can 
see that the formation of the liquid bridge is a very slow process, it is clearly diffusion-controlled, 
there is no evidence of faster hydrodynamic mechanisms. The radius of the midpoint of the bridge 
seems to grow towards its equilibrium value over a transient period of time according to a t 1 ^ 4 law. 
This behavior is reminiscent of the growth law with which an interfacial profile between coexisting 
phases approaches equilibrium 54 . 



5. Concluding discussion 

In the present work, we have presented Molecular Dynamics simulations addressing the shape of a 
droplet adsorbed on a spherical lyophilic domain on an otherwise lyophobic flat substrate surface, and 
the formation of liquid bridges between two such surfaces. Also the forces between these surfaces 
caused by such liquid bridges have been measured, and the kinetics of the merging of two such 
droplets into one bridge has been studied. 

According to the predictions of Lipowsky et al. 12 ~ 18 one should expect for large enough domain 
radius Rd that the contact angle 9 of the droplet is constant (9 — #p n jj, the "lyophilic" value, 
with cosflpj^j « 0.57 in our case, see inset of Fig. 2), until the radius R\ a ^ of the droplet matches 

Rd- For our choice of parameters, this happens for Rd = R D ^ ~ 12. For Rd < i?^'', one 
expects that 9 should increase such that R\ a i = Rd is always maintained, until 9 reaches the value 

of the lyophobic part of the substrate surface, 9 — #ph OD , for Rd = R^ ^ , with cosflpjjQb w —0.6 



in our case. In fact, Fig. 4 nicely verifies the linear variation i?j at = Rd quantitatively up to 
about Rd ~ 9, while for 9 < Rd < -R^^ the further increase of Ry^ with Rd is slower, and 

the saturation value of R\ & i, reached for R > iZ^"\ is in fact smaller than expected, namely only 
around ^™^ X ~ 10. Of course, some quantitative deviations of our results from the predictions of 
Lipowsky et al. 12 18 must be expected, since the latter predictions are asymptotically valid for very 
large, almost macroscopic, droplets, while our simulations concern nanodroplets, and hence there 
are no sharp transitions possible when Rd is varied: so Fig. 2 does not show a sharp kink of the 
curve cos(#) vs. Rd at Rd = R D ^ but rather a rounded crossover, and also at Rd = R^ ^ 

only a very smooth variation is seen (in fact, R^ ^ is of order unity in our case, and hence cannot 
be even uniquely identified from our data). Clearly, it would be interesting to analyze the amount 
of rounding of these "transitions" at 9 = #p n ;i and 9 — 0ph OD theoretically, and to repeat our 
simulations for larger droplets (and correspondingly chosen larger values of Rd) to see the extent 
to which these transitions become sharper, but all such extensions of our work would present major 
difficulties and hence have not been attempted. 

We also note that the theory implies that the shape of the droplet should be a sphere cap through- 
out. In particular, since the volume V of the droplet is constant, we should have 

^■(^lat + r2 ) = const 

and combining this equation with the geometrical relations for the contact angle , Eqs. (4), (5), 
in principle the variation of cos(#) with Rd in the regime R^ ^ 3 < Rd < Rjp 1 ^ can be explicitly 
predicted. Noting from Fig. 4 that for R D > R D ^ we have i?i t « 10, H « 7.5, one obtains 



V w 1400. However, for R D « 7 where 7? lat « 6.8, H w 9.3 (Fig. 4) we find that Eq. (6) would 
yield only V ~ 1100, and for Rd ~ 5, where R\ & ^ « 5, H m 11 we would get V ~ 1130, while for 
i?D = R\ a t = 3 where ss 13 we find V « 1330. Thus, the geometrical relations on the basis of the 
sphere cap picture are not verified accurately in our case. This problem is related to the fact that 

the linear increase of R\ & ^ with Rd does not continue all the way to R\ a ^ = R D ^ ~ 12 but R\ a ^ 
is significantly too small for Rd > 9. On the other hand, the discrepancies found are not dramatic 
either, they do not exceed 10-15%, in spite of the nanoscopic size of the droplets. One needs also to 
consider the fact that this small size leads also to substantial errors and ambiguities when one tries 
to read off R\at an d H from the data. In view of all these problems, the agreement between the 
simulations and the theoretical descriptions certainly is satisfactory. 
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FIG. 1. Contour diagrams representing the density profiles p(Z, R(X,Y)) = vAp, v = 1,2, ...11, 
Ap = 0.2, in the (Z, R) plane, for a droplet containing 128 chains with 32 monomers each, at a tem- 
perature ksT = 0.49. The droplet is in contact with an ideally flat lyophobic substrate (represented 
by the thick grey line in the bottom) decorated with one lyophilic circle of radius Rd (the thick 
black line in the bottom). The adsorption strengths of the lyophobic wall is ew = 0.05, while for the 
lyophilic circle it is ew = 0.2. Profiles are shown for Rd = 15, 12, 9, 7, 5 and 3, respectively. These 
profiles are obtained by averaging over 10 runs of 2.1 ■ 10 integration steps each. One MD time step 
is St = 0.0009 MD time units. 



FIG. 2. Cosine of the contact angle 9 plotted versus the radius Rd of the lyophilic domain, for 
droplets containing TV = 128 chains with N = 32 monomers each, at ksT — 0.49. The inset shows 
the dependence of the contact angle of the droplet at an infinitely extended homogeneous substrate 
(as considered in Refs. 37 39 ) on the strength ew of the wall potential. The two contour diagrams 
in the inset represent the density profiles p(Z, R(X, Y)) = vAp, v = 1, 2, . . . , 11 in the (Z, R) plane 
of a droplet for the two different adsorption strengths relevant for the present paper, namely for 
ew = 0.05 (the value used to model the lyophobic region of the wall) and for ew = 0.2 (the value 
used for the lyophilic domain on the surface). 



FIG. 3. Contour diagrams representing the density profiles p(Z, R(X, Y)) = fAp, v = 
1, 2, . . . , ll(Ap = 0.2) in the (Z,R) plane for two droplets adsorbed on opposite walls (L = 22) 
or a liquid bridge between the walls (L = 19, 17, 15, 13, 10, 9, 6, 5) respectively. The simulated sys- 
tem contains 256 chains with 32 monomers each, the temperature is chosen as ksT = 0.49, and the 
walls are ideally flat and lyophobic, with an adsorption strength e w = 0.05, except for a lyophilic 
circle of radius Rd = 14 (with adsorption strength e w = 0.2). The lyophilic parts of the adsorbing 
surface are shown by thick black bars, the lyophobic parts by thick grey bars. The profiles are 
averaged over 18 runs of 1.05 • 10 6 MD steps. 



FIG. 4. Variation of the base radius i?j at (circles), the height H (squares) and the midheight radius 
R\ & ^(Z = H/2) (diamonds) of a droplet adsorbed on a single circular lyophilic domain of radius 
Rd plotted vs. Rd at ksT — 0.49. The droplet contains 128 chains with 32 monomers each. The 
strength of the adsorption potential is e w = 0.2 inside the lyophilic domain and e w = 0.05 outside of 
it. The radii and the height are measured from density profile contour diagrams (taking the contour 
at midpoint density, p = 1) averaged over 10 runs of 2.1 • 10 6 integration steps each. R^ is taken 
as the maximum lateral extension of the contour (it occurs roughly at Z — 1.5). The error bars are 
obtained by estimating the radii from every single run. 



FIG. 5. Variation of the base radius i?j at (circles) and the midheight radius _Rj at (Z = L/2) (squares) 
for the liquid bridges of Fig. 3 (for further explanations see the caption of Fig. 4). A log- log plot of 
R^{Z — L/2) vs L (see inset) reveals a power-law relationship R\ & ^ cx L -0 ' 72 . 



FIG. 6. The pressure in a liquid bridge connecting two flat substrates is plotted against the distance 
L between the walls. The total pressure is shown by filled circles. Open diamonds show the contri- 
bution to the pressure that comes from the wall-monomer interactions, while the open squares show 
the contribution to the pressure from the monomer-monomer interaction and the kinetic part. 



FIG. 7. The normal force acting between two substrates with a bridging droplet plotted versus 
the distance L. The insets show contour diagrams of the density profiles of the liquid bridge at 4 
different distances L: the two separated droplets are at L = 22, the "in-bridge" is at L = 17, the 
cylindrical bridge is at L = 9, 5 and the "out-bridge" at L = 6. The thick dark lines indicate the 
lyophilic region of the substrate. Each point is obtained from averaging over 19 runs of 1.05 • 10 
MD steps, while the data points around the minimum (6.5 < L < 8.5) are averaged over 29 runs. 
The arrow at L = 22 indicates that bridges are no longer stable for L > 22. The contact angle at 
zero force is 9* « 100°. 



FIG. 8. Midheight radius R(Z — L/2, t) of a forming liquid bridge as a function of time, for L = 21.5. 
The radial distance from the vertical axis of the droplet to the point where the density is p — 1.0 



is measured. The origin of time t = is chosen as the time when the two droplets sitting on oppo- 
site parallel substrates have just touched each other. The insets show vertical cross sections of the 
resulting bridges for three different times indicated in the figure. R(Z = L/2,t) is averaged over 7 
independent runs. For 300 < t < 3000 an effective growth law R(Z = L/2,t) oc const + t ' 242 is 
observed, as indicated by a straight line. 



FIG. 9. Density profiles of a forming liquid bridge for four different times, as indicated in the figure. 
The diffusive interpenetration of species which originally belong to the two separated droplets is 
indicated. 
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